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1. Introduction 

The quantum integrable T2-model has been first introduced in [3 1 where its Lax operator was constructed as a 
general solution to the Yang-Baxter equations w.r.t. the 6- vertex R-matrix. There, for the subvariety of the cyclic 
representations parameterized by points on the algebraic curves associated to the chiral Potts (chP) modeU, the 
analysis of the spectrum (eigenvalues) was made by the standard construction of the Baxter Q-operator. This 
last operator was shown to coincide with the transfer matrix of the integrable chiral Potts (chP) model [4|-[8|; 
in this way a first remarkable connection between these two apparently very different models^ was established. 
The eigenvalues analysis was further developed in |[T2] [T3l where additional functional equations for the transfer 
matrices of these two models, like those corresponding to the fusion hierarch>§] of commuting transfer matrices^, 
were shown. In the special case of the superintegrable chP-model, the role of Bethe ansatz type equations for the 
spectral analysis was first pointed out in 15] |6]. In |[T6l the algebraic Bethe ansatz (ABA) description of the 
superintegrable chP-spectrum was rigorously introduced on the basis of the connection between the T2-model and 
the chP-model. In particular, the Bethe ansatz construction was applied to the T2-transfer matrix, which led to the 
reproduction of the Baxter results [81 on the subset of the translation-invariant eigenvectors of the superintegrable 
chP-model. More recently [17], Baxter has extended the eigenvalues analysis of the T2-model to completely 
general cyclic representations. The main tool used there was the construction of a generalized Q-operator which 
satisfies the Baxter equation with the T2-transfer matrixlf] and the consequent extension to these representations of 
the functional relations for the fused transfer matrices. 



1.1 Motivations 

Let us comment that, in the above literature, the spectral analysis suffers in general from at least one of the 
following problems: i) Analysis reduced to the eigenvalues only: no eigenstates construction for the functional 
methods based on the Baxter Q-operator and the fusion matrices, ii) Reduced applicability: the ABA applies 
only to very special representations of the T2-model and of the chP-model. iii) Lack of completeness proof, i.e. 
the completeness of the spectrum description is not given by the method. In particular, this last problem affects 
almost all the Bethe ansatz method^, which leads to the characterization of the spectrum by solutions of an 
associated system of Bethe ansatz type equations. 

A rigorous proof of the completeness of these spectrum characterizations is then a fundamental goal which was 
achieved in the literature only for a few examples of integrable quantum models, including for example the XXX 
Heisenberg model [ 32 ] and the cyclic representations of the lattice sine-Gordon model 11331 . It is worth remarking 
that in the case of the superintegrable chP-model and the T2 -representations to which the ABA applies such a proof 

'These curves are parametrized by the equations required for the periodicity of dilogarithm functions, i.e. (3.5) of 0. 

2 The differences between them are simply expressed in the 2-dimensional statistical mechanics formulation. Both are characterized by 
Boltzmann weights which satisfy the star-triangle (Yang-Baxter) equations but while those of the T2-model satisfy the ordinary difference 
property in the rapidities those of the chP-model do not. In fact, one of the interesting properties of the chP-model is that its weights 
are uniformized by curves of genus g > 1, hence they cannot satisfy the difference property. Let us recall that the first solutions of the 
star-triangle equations with this non-difference property were obtained in [9 10 1 1 1 , while in |7| the general solutions for the chP-model 
were derived. 

3 The approach of fusion hierarchy of commuting transfer matrices was first introduced in 111411151 . 

4 The T2-transfer matrix is the second element in this hierarchy, this explains the name given to the model. 

5 There, it was considered as a kind of generalized chiral Potts model transfer matrix, however, the fact that even the commutativity of 
the elements of such operator family is not given leaves the definition of generalized chP-model only formal. 

6 As the coordinate Bethe ansatz (I8][l9]|20), the algebraic Bethe ansatz (2T]-(29), the analytic Bethe ansatz I30II31I . 
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is missing. Instead, for the reduction to the trivial algebraic curve case, i.e. the Fateev-Zamolodchikov model 
11341 , the study of the completeness was addressed by a numerical method in ||35Tl (see also [36 ] and 11371 for 
further applications of this method). This stresses the importance of the proof of the completeness of the system 
of Bethe ansatz type equations for the T2-model and the chP-model that is given in this article. 

Note that all our results are derived in the framework of Skly anin's SOV method [38 ]-[40 ], which for the T2-model 
was first developed § in BT1 . However, we cannot get the spectrum characterization by functional equations 
directly from the SOV method. Indeed, in the case of cyclic representations [44] of integrable quantum models, 
this characterization is related to solutions of an associated finite system of Baxter-like equations. If this system is 
compatible with the eigenvalue equation associated to some Baxter Q-operator, the construction of this operator is 
one standard procedure to achieve a reformulation in terms of functional equations in the SOV characterization. 
For example, this was the strategy followed in ll33l for the special cyclic representations of the sine-Gordon 
model. In fiTTl . the functional relations of fusion and truncation of transfer matrices (generalized by Baxter 
in ifTTl ) were proven to be compatible with the SOV characterization, leading to the derivation of a functional 
equation reformulation for the T2 -model. 

Note that the standard Baxter construction of a Q-operator by gauge transformations and the subsequent deriva- 
tion of the fused transfer matrices can only be applied when the existence of some model dependent quantum 
dilogarithm functions P31 - ll55l is proven, which can represent a concrete technical problem. Then, it is relevant 
to ponder whether it is possible to bypass this kind of constructions by providing different derivations of the func- 
tional equation reformulation of the SOV characterization of the spectrum. We show that this is indeed possible 
for completely general cyclic representations of the T2-model; we achieve such a result following the approach 
first introduced in [1|. We prove that the T2 -transfer matrix (completed with the G -charge for some sub variety of 
representations) defines a complete set of commuting charges, i.e. it has simple spectruno Moreover, we provide 
the complete characterization of all its eigenvalues and eigenstates as solutions of given functional equations in 
well defined classes of functions. The same statements are proven for the inhomogeneous chP-transfer matrix. 
In the T2-self-adjoint representations, this characterization leads to the reconstruction of the T2-eigenbasis which, 
for the reduction to representations on the chP algebraic curves, are proven to be a simultaneous eigenbasis of the 
inhomogeneous chP-transfer matrix. Furthermore, we show that there exists a quite general subvariety of T2 -self- 
adjoint representations, for which we can reconstruct the eigenvalues and the full basis of eigenstates in terms of 
the solutions of an associated system of Bethe ansatz type equations. Finally, we prove the same statement of 
completeness for the chP-spectrum, which corresponds to the special T2-self-adjoint representations on the chP 
algebraic curves. 

In the paper |[82l an approach has been developed in the framework of the quantum inverse scattering method 
(QISM) 112111 - 11291 , [83]-] 86] to achieve the complete solution of lattice integrable quantum models by the exact 
characterization of their spectrum and the computation of the matrix elements of local operators. In particular, 
in [|82 ] the approach has been developed for the relevant case of the lattice quantum sine-Gordon model 11241 l84l 
associated by QISM to some special cyclic representations [44] of the 6- vertex Yang-Baxter algebra, and it can 
be considered as the generalization to this SOV framework of the Lyon group method^. The possibility to adapt 

7 See also the series of works 1142 1 1431 where a first result concerning the computation of the form factors of local operators by SOV 
was achieved for the special case of the generalized Ising model. 

8 It is worth remarking that the simplicity implies that the fused transfer matrices do not really add further information useful for the 
spectrum characterization. 

9 Let us remark that this method has been introduced in [87 1 in the algebraic Bethe ansatz (ABA) framework 1211 - 1241 for the spin-1/2 
XXZ quantum chain 1 18 1, [88 1-[95 1 with periodic boundary conditions and then further developed in 1961 - 111081 . See also 1 1 09 1 1 1 1 01 and 
11111 - 11131 . for the extension of this method in the ABA framework, to the higher spin XXX quantum chain and to the open spin-1/2 
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this approach to the T2-model using the SOV framework is another motivation to this paper. 



1.2 Organization of the paper 

The paper is organized as follows. In section 2 we introduce the T2-model for general cyclic representations, 
pointing out its main properties and characterizing the representations for which its transfer matrix is self-adjoint. 
In section 3 we characterize the spectrum (eigenvalues/eigenstates) for the most general representations by using 
the quantum separation of variables only. In particular, we prove the simplicity of the T2-spectrum and we show 
how to reformulate the SOV characterization of the transfer matrix spectrum in terms of solutions of a functional 
equation. In section 4 we consider the restriction of the T2-model to self-adjoint representations for which we 
show that a propeiF"! Baxter Q-operator is naturally induced from our SOV-characterization of the T2-spectrum. 
Moreover, we characterize the most general T2 -representations for which such a Baxter Q-operator is polynomial 
in the spectral variable, leading to the proof of the completeness of the associate Bethe ansatz equations. In 
section 5 we introduce an inhomogeneous version of the chiral Potts model and we characterize all its eigenstates 
by those constructed for the T2 -transfer matrix. Moreover, we prove the normality of the chiral Potts transfer 
matrix for the self-adjoint T2 -representations on the chP algebraic curves, which completely characterizes the 
chiral Potts spectrum (eigenvalues/eigenstates) in terms of the one of the T2-model. Finally, we prove that there 
exists a nontrivial subvariety of these representations for which the simultaneous spectrum of chiral Potts model 
and T2 -model is completely characterized by polynomial solutions to the associated Baxter equation also proving 
a completeness statements for the solutions to Bethe ansatz equations. In appendix A, we give, in some details, the 
construction of the SOV representation for the T2-model. In appendix B, we reproduce, adapting to our notation, 
the Baxter construction of the generalized Baxter Q-operator for the T2-model, to point out the connections with 
our SOV construction of the ti -spectrum. Finally, appendix C describes in the most general framework of cyclic 
representation some useful technical results. 

2. The r 2 model 

2.1 Definitions and first properties 

The Lax operator which describes the T2 -model can be parametrized as followsB 

L [X)= { \ot n y n - PnX- 1 ^- 1 u„ (<T 1/2 a„v n + q 1/2 h n y- 1 ) \ ^ ^ 

\ u,; 1 (g 1/2 c n v n + g _1/2 d n v~ 1 ) 7nV n /A - £„A/v„ y' 

where the a„, h n , <c n , <d n , a n , (3 n , j n and 6 n are constants associated to the site n that satisfy the relations: 

a„7„ = a„<c„, f3 n S n = b„d„, (2.2) 

and we have denoted with u„ and v„ the generators of the local Weyl algebras: 

u„v m = </" m v m u„ Vn,m £ {1, N}. (2.3) 

We will study the case where q = e ™ p is a p-root of unity: 

/3 2 = - , p = 21 + l,p' = 21' and I, I' G Z >Q -> <f = 1. (2.4) 
V 



XXZ quantum chain [ 1 14|-| 120 1 with diagonal boundary conditions, respectively. 

10 Here we are pointing out that this Q-operator define a one-parameter family of commuting operators which commute with the T2- 
transfer matrix. 

1 'Up to different notations, this Lax operator coincides with that introduced in 1 3 1. 
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In this case we can define a finite-dimensional representation of dimension p for each Weyl algebra W n by introducing the 
states: 

\z) = \ Zl ,...,z N ) with z t G Bp = {q 2n ;n = 0, ...,21} and i e {1, N}, (2.5) 

for which we have: 

v„ \zi, ...,zu) = I zi, . . .,qz n , ...,z N ), 
u„ \zi,...,z N ) = z n \ zi,. ..,z N ). 
From the Lax operators, we can define the monodromy matrix by: 



(2.6) 



which satisfies the quadratic relation known as the Yang-Baxter relation: 

R{\/p) (M(A) ® 1) (1 ® M(/i)) = (1 ® M(/i)) (M(A) ® l)-R(A//x) , (2.8) 
w.r.t. the six-vertex (standard) i?-matrix: 



V qX-q^X- 1 ) 



(2.9) 



The elements of M(A) then generate a representation 7\Ln of dimension p N of the so-called Yang-Baxter algebra, this 
representation is characterized by 6N parameters. In particular, the commutation relations (12.81 l lead to the relation 
[B(A), B(/i)] = for all A and \i, and also to the mutual commutativity of the elements of the one-parameter family of 
operators: 

t 2 (A) ee tr c2 M(A) = A(A) + D(A), (2.10) 
known as transfer matrix. Let us introduce the operator: 

N 

<d=Y[v n , (2.11) 

n=l 

which plays the role of a grading operator in the Yang-Baxter algebrcQ 

Lemma 1. commutes with the transfer matrix. Moreover, it satisfies the following commutation relations with the ele- 
ments of the monodromy matrix: 

6C(A) = gC(A)0, [A(A),6]=0, (2.12) 
B(A)0 = g0B(A), [D(A),0]=O. (2.13) 

Besides, the O-charge allows to express the asymptotics of the transfer matrix in A — > and in A — > oo; in particular, from 
the known form of the Lax operator, we derive the following expansions: 

/ N N \ N-l 



A(A) = A N 9 J] a n + (— 1) N A~ N 1 I] ^ + E A ' AN_2t ' i2A ^ 

\ a—1 a—1 / i—l 

/ N N \ N-l 

d(a) = A- N e [] 7a + (-i) N A N e- 1 [] s a )+ £ d * aN ^ (2- 15 > 



a=l 



with Aj and being operators, and s 



lim A ±N r 2 (A) = (6 Tl a T + ±1 d T ) , (2.16) 

log A— ^=poo 



12 The proof of the lemma is given following the same steps of that of Proposition 6 of [331. 

13 Here, we have used the short notation limiog a->±oo, where — stands for the limA->o and + stands for the lim> 
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where: 

N N N N 

a+ =n aQ , a _^(-i) N n^ a , d+=(-i) N jjd a , d_=jj 7a . (2.i7) 

a— 1 a— 1 a— 1 a— 1 

Let us denote by £ T2 the set of the eigenvalues t(X) of the transfer matrix T2(A), then: 

S T2 C C CTen [A, A _1 ] N for N even, E T2 c C odd [X, A _1 ] N for N odd, (2.18) 

where C e [x, x _1 ]m is the linear space of the Laurent polynomials over the field C of degree M in the variable x that are 
even or odd, as stated in the index e. The O-charge allows to introduce the grading £ T2 = UfcLo ^r 2 ' wnere: 

S^ 2 ee (t(A) e S T2 : lim \~ N t(\) = (q^ k a T + q ±k d T )\ , (2.19) 

L logA^=poo J 

where to any t(A) £ Sj: 2 corresponds simultaneous eigenstates of T2(A) and with 0-eigenvalue 



2.i.i Quantum determinant 
The quantum determinant: 

det q M(A) ee A(A)D(g _1 A) - B(A)C(g- 1 A), (2.20) 
is a central element of the Yang-Baxter algebra ( 12. 8t which has the factorized form: 

N 

detqM(A) = JJ detqL„(A), (2.21) 
n=l 

where: 

det q L„(A) ee (L„(A)) n (L n (A/g)) 22 — C-n)i 2 (L n ) 2 i , (2.22) 
are the local quantum determinants, which explicitly read: 

det q M(A) = \{k n (^- - - ^) 

n— 1 ' ' 

= (-9) II (t+9 — B-A)(-+ g — —A), (2.23) 

„ i a n A a„p„ A <n n p„ 

n— 1 

where: 

j / a i a/2 J «9 1/2 (a„/3„/a n b n ) 1/2 ft, = +, 

fe„ = (a„b„<E„<d n ) , = < . , 1/2 (2.24) 

[ [€ n p n / a n d n ) h = -. 

2.2 Self-adjoint representations 

Lemma 2. Lef e G {—1,-1-1}, if the parameters of the representation satisfy the constrains: 

c„ = -eb*, d„ = -ea* , /3 n = e (a*b n ) /a* , (2.25) 
then the generators of the Yang-Baxter algebra obey the following Hermitian conjugation relations: 

M(A) = \0{X) Dt(A) J = (,-eB(A*) A(A*) J' (2 ' 26) 

which, in particular, imply that the transfer matrix t 2 (A) is self-adjoint for real A. Besides, the quantum determinant has 
the expression^\ 

, , y , u N TT |a„| 2 |b„| 2 1 -i|an| 2 ,wl . _i K| 2 u , ™ 

det q M(A) = , JJ___(- + e « KF A)( A +£,? jto^TF ( } 



14 See 1 56 1 and also 1571 for an historical note on the centrality of the quantum determinant in the Yang-Baxter algebra. 
15 Remark that it only depends on the modulus of the parameters of Lax operators. 
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Proof. It is simple to observe that the Lax operator of the -model satisfies the equation ( 12.26b . which can be also written 



as: 



(l n (\)) f = <r 1+5l , e L„(A*) a 1+Sl c — > M(A)t = a 1+Slt M(A*) a 1+Sl c , (2.28) 
that is d2T26l > holds by definition of M (A). □ 

3. Characterization of r 2 -spectrum: general representations 
3.1 SOV representations 

According to Sklyanin's method [38 39] S3, a separation of variables (SOV) representation for the spectral problem of 
T2(A) is given by a representation where the commutative family of operators B(A) is diagonal. 

Theorem 1. For almost all the values of the 6N parameters of the representation, there exists a SOV representation for the 
T%-model, i.e. B(A) is diagonalizable and has simple spectrum. 

Proof. See appendix|A]for a constructive proof of this statement. □ 
Let ( 77k | be the generic element of a basis of eigenvectors of B(A): 

N-l 

( Vk |B(A) = 6, k (A) (r /k | , ^(A) = [] (V^f * 3 ~ , (3-D 

0=1 

and 

ry k G Z B = { ( V [ kl) ^ q kl V ( i\- ■ • , = q""^) > k = ■ • ■ . M * K } . 0-2) 

where rf^ are constants that are defined in appendix [A] Here, the simplicity of the spectrum of B(A) is equivalent to the 
requirement {jlaj ^ (jlb*^ f° r an y a^bg {1,...,N- 1}. The action of the remaining generators of the Yang-Baxter 
algebra on arbitrary states^ ( 77 1 then reads: 

where af- sov ^ (ry a ) and d( soy ' (rj a ) are coefficients which have to satisfy the quantum determinant condition: 

det q M(f>) = a (soy) (v>)d (soy) (^V) , Vr = 1, . . . , N — 1 . (3.5) 

Here, we have defined: 



N-l N-l 

Va 



(±) = (±i) N - i o ± n nt\ t& = (±i) n -^ ± n ^ (3.6) 



n— 1 n— 1 



and the states ( q ± a rj | are defined by: 

(q ±s °r l \ = (r ]1 ,...,q ±1 r ]a ,..., m \. (3.7) 
Finally, C(A) is uniquel}0 defined by the quantum determinant relation. 



16 From here on to simplify the notation, we will omit the subscript k in 77k as well as the superscript (k a ) in the rja and we will 
reintroduce them only when it will be strictly required. 

17 Note that the operator B(A) is invertible except for A which coincides with a zero of B, so in general C(A) is defined by (4.5) just 
inverting B(A). This is enough to fix in an unique way the operator C as it is a Laurent polynomial of degree N — 1 in A. 
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3.2 Average values of Yang Baxter generators as central elements 

We define the average value O of the elements of the monodromy matrix M(A) by: 

0(A) = J] 0(q k X) , A = X p , (3.8) 

k=l 

where can be A, B, C or D and the commutativity of the families A(A), B(A), C(A) and D(A) implies that A(A), T>(A) 
are Laurent polynomials of degree N while B(A), C(A) are Laurent polynomials of degree N — 1 in A. Then the following 
proposition holds: 

Proposition 1. 

a) The average values „4(A), B(A), C(A), T>(A) of the monodromy matrix elements are central. Besides, they satisfy the 
relation: 

(.4(A))* = 2?(A*), (S(A))* = -eC(A*), (3.9) 
in the case of self-adjoint representations. 

b) LeB 

^ = ( $ r, ) 

be the 2x2 matrix whose elements are the average values of the elements of the monodromy matrix M(A), it holds 
then: 

M(A) = i C N (A)£ N _ 1 (A) ...A(A), (3.11) 

where: 

Ln{A) -[ q ^« + o Tg/A-A^: ) (3 ' 12) 

is the 2x2 matrix whose elements are the average values of the elements of the Lax matrix L n {\). 
Proof of a). 6(A) is central as it follows by taking directly the average of (13.11 ) in the SOV representations: 

N-l 

B(A) = Z N Y[ (A/Z a - Z a /A) , Z a = v P= [rji^Y . (3.13) 



From q p = 1, we have that A(Z r ) and T>{Z r ) are central and related to the coefficients aS sov > (q k r] r ) and d^ sov \q k rj r ) 
by 

A(Z r ) = f[^ so ^(q k n r ), V{Z r ) = f[a( sov \q k Vr), Vre{l,...,N-l}. (3.14) 

k=l k=l 

^l(A)A N and 2?(A)A N are polynomials in A 2 of degree N. The relations ( 13.141 and the simplicity of the B-spectrum give N 
points in which these polynomials are central elements, and the centrality of the asymptotics of .4(A) and P(A), as trivially 
follows from (l2TT4l - ([2T3t . yields the centrality of .4(A) and V(A). □ 

Proof of b). By using that B(A) is diagonalizable and with simple spectrum in the entire chain as well as in each subchain, 
the point b) follows inductively by the simple extension to our representations of the recursion relations on the averages of 
the Yang-Baxter generators of Proposition 3 of ll33l . 

□ 



8 A similar statement was first proven in 11441 . 
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3.3 SOV characterization of r 2 -model 

The spectral problem for T2(A) in the SOV representations is equivalent to the discrete system of Baxter-like equations for 
the wave-function ^t(v) = ( V 1 1 ) °f a 72-eigenstate 1 1 ): 

t(t] r )Mv) = ^ sov) M^t(q- Sr v) + <i {SOV) (Vr)'ft(q Sr v) Vr6{l,..,N-l}, (3.15) 
plus the equation in the variable rjfi : 

Mq*"v) = q~ k vM, (3.16) 

fort(A) e with k e {0, ...,21} where (r?i, r] N ) € Z B and 

q ±5r j] = (t]i, . . . ■ ■ ■ ,vn)- (3.17) 

Let us remark that the coefficients a^ ^ (rj r ) and dS sov ' (r) r ) which satisfy d3.5t and the average conditions ( 13.14| > are 
fixed up to the following gauge transformations: 



(3.18) 



f(Vr) ' " V "' " V '" /(»7r) 

which just amounts to a renormalization in the states of the B-eigenbasis: 

N-l 

Hr'Mivl- (3.19) 

r=l 

Here, we make the following choice of gauge: 

a (S ° V) (A) = ( TTF^LrV /P A(A), d^>(A) = ( rrp^^ Y'vW, (3.20) 

where: 



i1q=i A (M a )j w ' w ~vnLiD(A? a 



N . N 



a (A) ee TT(^ an ) 1/2 (^- - d(A) ee TT( a " b "^ d " )i/2 _ t^y (3.21) 

— 1 n— 1 



Remark 1. The average values Z r , Zj^, Zq, A{Z r ) and V(Z r ) are cleary unchanged by gauge transformations. Moreover, 
as central elements of the representation, they are not modified by similarity transformations, and they characterize gauge- 
invariant parameters of the SOV representation. Moreover, Proposition [T] uniquely defines the average values of the mon- 
odromy matrix elements and allows to establish that .4(A), £>(A), C(A) and T>(A) are polynomials of maximal degree 1 in 
each of the parameters a v n , (3?, a v n , 6^, c^, d v n of the t-i -representation. Note that this also implies that the gauge-invariant 
dates of the SOV representations (up to permutations) are uniquely defined in terms of these parameters. In particular, 
denoting by a„ 1 (Z) the elementary symmetric polynomial of degree n in the var iables Z r , the a^' 1 (Z) / a^Z\{Z) are 
polynomials of degree 1 in the parameters of the T2-representation. 



3.4 Simplicity of T2-spectrum 

In this section we show that the spectrum of the transfer matrix (A) is non-degenerate (or simple). Let us prove that: 

Lemma 3. For almost all the values of the parameters (pP n , f3f ll afj, fo^, c^, d£) the average values of the monodromy 
matrix elements A(A) and D(A) satisfy the inequalities: 

A(Z a )^V{Z a ), Vae{l,...,N-l}, (3.22) 

where Z a are the zeros of the average value of B(A). 
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Proof. Let us define the functions^ 

T a {al, (3P, a£, K, ^, d^) = A N (Z a ) - V N (Z a ), Va G {1, N - 1}. (3.23) 

The Proposition Q] and the Remark 1 about the functional dependence of Z±, . . . , Z N _ X w.r.t. these parameters implies that 
it is sufficient to show that the functions T a are nonzero for some special value of the parameters in order to prove that they 
are nonzero for almost all the values of the parameters , [3?, b£, c^, d£, which will prove the lemma. Note that the 
following identities hold: 

„ Ii9l = wi* ( C? + <)^K + K) =g p/2 ((^ + 0-(a? + b?))B 2 , N _ 1 (Z ), (3.24) 

( d„ = a„, E n = b„, /3„ = -b„a„/a„, n£{2 N}, 

Here, we have used the decomposition of the chain in a first subchain l, formed by the site 1, and a second subchain 2, 
formed by the remaining sites. We also have used the identities: 

(C ltl -B ltl ) = <f/ 2 (« + d?)-( a ? + b?)), A,i(A)^u(A), (3.25) 

B 2) n-i(A) = C 2 ,n-i(A), D 2)N _ 1 (A) = A,n-i(A), (3.26) 

which follows by using Proposition Q] for the special choice of the parameters done in ( 13.241 >. Then, we have only to show 
that we can always take B Z: n-i(Z a ) ^ for any a G {1, N — 1}. 

Note that from the identities: 

det q X 2 , N _!(A) = (A,n-i(A)) 2 -(S 2 ,n-i(A)) 2 , (3.27) 
B N (A) = A,i(A)S 2 , n _ 1 (A) + 6 1j1 A,n-i(A), (3.28) 

we have that B z , N-i(Z a ) = if and only if Z a is a double zero of det q A^ 2i n-i(A). However, this is not the case in 
general: by averaging the quantum determinant (12.23) , we get the formula: 

N 

det q X 2 , N _i(A) = J] <K II ( A / M n,h - M n , fc /A), M n ,+ = a£/<, M„,_ = b£/<, (3.29) 

so the function det q .M 2j n_i (A) has not double zeros for general values of the parameters a„, b„, a„. □ 

Theorem 2. For almost all the values of the parameters a£, a^, b£, c^, d^ of a ^-representation, the spectrum of 
ra(A) is simple. 

Proof. We have to prove that, up to normalization, for any given t(X) G S T2 there exists only one solution to the system 
(13 -15b and ( 13.161 1. Let us denote by ^t(v) w ith rj G Zb any other solution corresponding to the same r 2 -eigenvalue t(X). 
Then, we can define the q-Wronskian: 

WtAv) = Mv) *t(<rS) - Mv) *t(«~N), vr g {i, n - 1}, (3.30) 

which by the system of Baxter-like equations (13.151 satisfies the equations: 

^ sov \vr)W t Av) = d (SOV) (Vr)W t>r (q Sr V), Wg{1,...,N-1}. (3.31) 
Thanks to the cyclicity the averages of the above equations read: 

22 

(A(Z r )-V(Z r ))WtAv) =0 with WtAv) = Il^rfe 9** In), (3-32) 

fc=0 

which by Lemma|3]implies: 

W t Av) = 0, V?; G Z B and Vr G {1, N - 1}. (3.33) 



"From here, we will use the index N when it will be needed to point out that we are referring to the chain with N sites and we will 
omit it otherwise. 
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Note that this implies that the wave functions are proportional to the following factorization: 

N-l 



(3.34) 



where the proportionality factor could depend on fixed parameters, such as ^i "* j , and the Qt(i]r) is the solution of the 
discrete system of Baxter-like equations (13.15b for the fixed r e {1, N — 1}. □ 



3.5 Characterization of T2 -eigenvalues as solutions of a functional equation 

Let us define the following functions: 

a(A) ee a(A)A(A), 13(A) = ar x (gA)D(A), (3.35) 
where a(A) and d(A) are the Laurent polynomials defined in (13.211 1. Then, they always satisfy the condition: 

det q M(A) = a(A)d(A/<?), (3.36) 
while the function a(X) is defined by the requirement: 



Yl A(Xq n ) + f[ 5(A? n ) = A(A) + 23(A). 



(3.37) 



Note that this last condition is a second order equation in the average n«=i o>(q n X) and then we have only two possible 
choices for the averages of the functions a(A) and 5(A): 



II A(\q n ) = il e (A) , n D(Ag") = n_ e (A) , 



(3.38) 



n=l 



where e = =F and Vt± are the two eigenvalues of the 2x2 matrix A4 (A) composed by the averages of the Yang-Baxter 
generators. Let us introduce the one-parameter family D(X) of p x p matrix: 



( t(A) -5(A) 
-A(gA) t{q\) -D(flfA) 





D(X) 







-A(A) \ 






\~D(q 2l X) 





A^-U) t{q 2l ~ l X) -f>{q n - l \) 
~A(q 2l X) t{q 2l X) J 







(3.39) 



where for now t(A) is only a Laurent polynomial of degree N in A, even for N even and odd for N odd. 

Lemma 4. The determinant of the matrix D(X) is a Laurent polynomial of maximal degree N in A ee X p , even for N even 
and odd for N odd. 

Proof. Let us start observing that D(Xq) is obtained from D(X) by exchanging the first and p-th column and then the first 
and p-th row, so that 

det D(Xq) = det D(X) VA e C, (3.40) 

p p 

then detp D(X) is a function of A. Let us observe now that dct p 13(A) admits the following expansion: 
detD(A) = -(^(A)+23(A))-A(A)D(A/g)detD (1>2m3>(1 , 2m) (A) 



21-1 



-A(Ag)D(A) det r>(i,2),(i 1 3)(A) + t(A)detDi ll (A) P 



(3.41) 
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where -D^^w^m (A) denotes the (21 — 1) x (21 — 1) sub-matrix of D(X) obtained removing the rows and columns h and k. 
The tridiagonality of the matrices Di,i, -0(1,2), (1,2)) ^ (i,2i+i),(x,2j+i) implies that their determinants coincide with those 
of the matrices obtained substituting the functions A(A), d(A) with the Laurent polynomials a(A), d(A). Then, the lemma 
follows as all the terms in the expansions ( 13.4 U are Laurent polynomials with the same properties stated in the lemma. □ 

The interest toward the function det p D(A) comes from the following result: 

Lemma 5. Anyt(X) G S T2 is a solution of the functional equation: 

detL>(A)=0. (3.42) 
p 

Proof. First note that the determinant det p D(A) depends from the coefficients A(/x) and E>(/i/g) only through their products 
(13.361 1 computed in fi = q a X with a = 1, ...,p (i.e. by the quantum determinant) and the sum of their averages (13.37b 
computed in A (i.e. by the sum of the averages of the operators A and D). Indeed, this statement trivially follows from the 
expansion (13. 4U and from the tridiagonality of the matrices Di,i, -On 2), (1,2)) -Dji^i+i), (1.22+1)- 

From the previous lemma, we only have to show that the determinant is zero in N — 1 points, and that the asymptotics in 
±00 are also zero. Let us observe that the SOV characterization of the T 2 -spectrum implies that the system of equations 
(13.15b admits a non-zero solution, i.e. t(X) G £ T2 only if: 

detD(r£)=0 Va G {1, N — 1} and (r?i, r) N ) G Z B . (3.43) 
p 

Indeed, the SOV coefficients lead to the same values of the quantum determinant and the averages of A and D in the zeros 
of B. Besides, we have: 

lim A ±N detD(A) = 0, (3.44) 

log A— >=foo p 

which simply follows by observing that: 

lim A ±N detL>(A) = - det \\a=5i ,_i + d=6 ij+ i - (q Tk a T + q^d^S.Al = 0, (3.45) 

log A— >=poo p p " ' ' " 

fori(A) G Y, k T2 andfc G {0, ...,21}. □ 

Remark 2. The same type of functional equation dct£>(A) = also appears for different quantum integrable models 
in ll58l l59l l60l . There, it stands for the functional relations which result from the truncated fusions of transfer matrix 
eigenvalues. In the case of the T2-model this type of fusion leads to the same type of equation, as it has been derived in 

|H2 [12 El El 

3.6 Construction of T2-eigenstates 

Thanks to the previous results we can give a complete characterization of the set E T2 and construct one r 2 -eigenstate \t) for 
any t(X) G S r2 . 

Theorem 3. S T2 coincides with the set of solutions to: 

detL>(A)=0, VAgC, (3.46) 
p 

in the class of functions ( I2.7<SI ). Then, we can associate a T2-eigenstate: 

N-l 

* t fa) = (m,-,VN \t) = v N k n 0.47) 

r=l 

to any t(X) G S^ 2 , where Q t (A) is the unique solution ( up to quasi-constants) corresponding to t(X) of the Baxter equation: 

t(X)Q t (X) = A(X)Q t (X/q) + D(X)Q t (qX). (3.48) 
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Proof. Lemma implies that any element of E r2 is a solution of (13.461 1. The reverse inclusion is given by the fact that the 
condition (|3.46l l is just the requirement of the existence of a non-zero solution Qt(X) of the equation ( 13.481 ). Then, t(X) 
and the ^t(v) defined in (|3.471 i are solutions of the discrete system of Baxter-like equations d3.15t and d3.16t and so they 
define a T2-eigenvalue and the corresponding T2-eigenstate. Let us show that Qt(X) is unique; denoting with Qt(X) any 
other solution, we can define the q-Wronskian: 

W t (X) = QtWQtiq-'X) - Qt(A)g t (q- 1 A) , (3.49) 

which by the Baxter equation satisfies the equation: 

A(X)W t (X) = b(X)W t (qX). (3.50) 

Thanks to the cyclicity the average of the above equation reads: 

P P p 

(l[A(Xq k )-l[D(Xq k ))W t (A) = with W t (A) = JJ W t (Xq k ), (3.51) 

fc=l k=l fc=l 

which implies W(X) = 0, being: 

p p 

Yl A(Xq k ) + Yl 5 (V°) or equivalently fl+ (A) ^ (A) . (3.52) 

fc=i fe=i 

It is then easy to see that this implies that Qt(X) = Qt(X) up to quasi-constant normalization. □ 

4. Characterization of r 2 -spectrum: self-adjoint representations 

4.1 SOV reconstruction of a Baxter Q-operator 

The interest toward the self-adjoint representations of the T2-model is due to the use of the spectral theorem. The ra(A) 
transfer matrix is then diagonalizable and the Theorem [3] leads to the construction of the T2-eigenbasis; this in particular 
means that the self-adjointness of T2(A) forces the Baxter equation (13.48b to have a complete set of independent solutions. 
Moreover, it is interesting to note that in the self-adjoint representations, from the results of our SOV analysis, a Baxter 
Q-operator for the T2-model is automatically reconstructed. In particular, we can define the operator: 

Q(X)\t) ~ Q t (X)\t) (4.1) 

where Qt(X) satisfies the Baxter equation (13.48b with t(A) £ E t and |i) is the corresponding T2-eigenstate. Then, the set of 
the T2-eigenstates being complete and the function Qt(X) being uniquely defined (up to quasi constants), the operator Q(A) 
is well defined, up to a quasi-constant scalar quantity, and clearly satisfies the properties: 

Q(A)r 2 (A) = A(A)Q(A/g) + D(A)Q(Ag), [Q(A),Q(/x)] = 0, [r 2 (A), Q(/i)] = VA^eC. (4.2) 

In the next subsections we restrict our attention to special self-adjoint representations of the T2 -model. Similarly to the case 
of the sine-Gordon model described in 11331 Q] |2], we show that for these self-adjoint representations the transfer matrix 
spectrum (eigenvalues and eigenstates) is completely characterized in terms of polynomial solutions of the associated 
functional Baxter equation. Moreover, we prove the completeness of the set of the transfer matrix eigenstates constructed 
by using SOV from the solutions of the associated Bethe ansatz equations. 

4.2 Construction of polynomial Baxter equation solutions from t 2 -eigenvalues 

Let us consider the subvariety in the self-adjoint T2 -representations of real dimension 4N characterized by the constrains: 

fj< =h ^ = ^ ; Vne {!,..., N}. (4.3) 

a h b* a* a n+ ia n b„+i&* 
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Lemma 6. The Laurent polynomials: 

a(A) ee ^ N f[ %(1 " i^/V^H^ 1 - i^/V 1 / 2 ^), d(A) ^ g N &( _ Ag)) (4 4) 

x \ A a„ b„ 

where e = ±, satisfy the equations: 

det ? M(A) = a(A)d(A/g), VA e C, (4.5) 

one/ 

p p 

.4(A) + 2?(A) = Y[ a(Ag n ) + JJ d(Aa n ), VA e C. (4.6) 

n— 1 n— 1 

Proof. The condition ( 14.5b is trivially verified. The condition (|4.6l l is instead a consequence of Proposition [6] and, in 
particular, of the average properties (|B.63) - (lB.64t once we prove that the coefficients a(A) and (Zb(A), defined in ( IB.38I ), 
have the same average values. This last statement trivially follows by direct verification once we make the identifications: 



b\l/2 / jjr \l/2 

n-ia«-i \ _ /_ a n a, n N 

K_l<_l/ V a " a n 



Vne{l,...,N}, (4.7) 



which are consistent with the requirement of cyclicity (IB. 22b under the constrains (14.31 l. □ 

Remark 3. Let us remark that the lemma can also be proven under the condition: 

A(A)=V(A), (4.8) 



and that the self-adjoint T2-representations for which (14.8b holds also in any quantum site, i.e. those for which 
(ai, qn, ai, ai\i, hi, bpj) S R m , is a subvariety of the one defined in (14.3b . 

It is worth remarking that, in the class of self-adjoint representations defined by the constrain ( 14.3b . it results: 

a + = d + and a_ = cL_. (4.9) 

Then, in these representations, the transfer matrix spectrum presents a double degeneracy for |fc| ^ 0, since the asymptotic 
behaviour of a i(A) € 2-v 2 and a i'(A) £ coincide. Such a degeneracy in the T2-spectrum is resolved by the 0-charge, 
i.e. the couple (t2(A), 9) has simple spectrum. The proof of this last statement can be given by adapting the proofs of 
Proposition 5 and Lemma 2 of f33 1 to these self-adjoint representations of the T2-model. 



Theorem 4. Let us take a t(A) G S^ 2 , then there exists, up to multiplicative quasi- constants, a unique polynomial of the 
form: 

2iN-(bt+o«) 

Qt(A) = A Q( J] ( A ^ A )- °<«t < 2 ^ 0<6 4 + at <2/N, (4.10) 
which satisfy the Baxter equation: 

t(X)Q t (X) = a(A)Qi(Ag _1 ) + d(\)Q t (Xq) VA 6 C, (4.11) 

vv/zere: 

at = ±fc modp, ot = ±fc mod p. (4.12) 

Proof. The uniqueness up to multiplicative quasi-constants is a consequence of Theorem[3] So we only have to prove the 
existence of such a polynomial Qt(X) and the fact that it has the form defined in (14.10b and (14.12b . First, let us define the 
matrix D(X) by replacing A and D by a and d in D(X). From the Lemma|6]and since the determinant of D(X) only depends 
of A and D via their average values and the quantum determinant, D(X) and D(X) have the same determinant. It is worth 
noticing that the condition t(A) G Sj: Q implies that the p x p matrix D(X) has rank 21 for any A e C\{0}. Then denoting 
with 

Ci d {X) = (-iy +1 det (4.13) 



17 



the cofactor of the matrix D(X), the matrix of elements these cofactors has rank 1 and so the vectors: 

Vi(A) = (C M (A), C il2 (A), Q, 2/+1 (A)) T eP Vi e {1, 2/ + 1} (4.14) 
satisfy the proportionality conditions: 

V i (A)/C M (A) = V i (A)/C J - 1 (A) Vi,je{l,...,2Z + l}, VA g C. (4.15) 

These conditions implies: 

C 2 ,2(A)/C 2il (A) = Ci,a(A)/Ci,i(A), (4.16) 

and for (IC.lt we can rewrite it as: 

Ci J i(A5)/Ci, a+ i(A«) = C 1)2 (A)/C 1)1 (A). (4.17) 

Moreover, from the relation: 

Ci,2J+i(A) = q N C h2 (-X/q), (4.18) 
which follows from dC.lt and (IC.2K (14.16t is also equivalent to: 

Ci,i(Ag)Ci,i(A) = g N C 1 , 2 (A)C 1 , 2 (-A) = gn N Ci,2i + i(A)Ci,2i+i(-A), (4.19) 

and from the condition det p D(A) — 0, we have: 

*(A)Ci,i(A) = a(A)Ci, M+ i(A) + d(A)Ci, 2 (A). (4.20) 

Let us note that all the cofactors are Laurent polynomial of maximal degre j@2ZN in A: 

4/N-(<H, 3 +6i,j) 

C itj (\) = C id \- 2m+a ^ \[ ' (AP-A), (4.21) 

h=l 

as it follows from the form of a(A), d(A) and i(A) G S T2 . Thanks to formula (IC.2t . the cofactor Ci i(A) G C[A, A _1 ] 2 ;n is 
even in A: 

2iN-(5i,i+6i,i) 

C 1 , 1 (A) = c 1 , 1 A- 2iN + 25 ^ [] (Af' 1) -A)(A| 14) +A), (4.22) 
then the cofactor equation (14.19t implies: 

ai,2 = 2oi,i = 2a, 6i, 2 = 26 M = 26, c 2 2 = c 2 hl q- 2{N+h} (4.23) 

and: 

(Af^) 2 =(AP) 2 .Af, (AS^) 2 ^/,) 2 , (4-24) 
with Ai 7^ for any i G {1, 2ZN — (a + b)} with a and b G Z-°. Then, we can write: 

2/N-(a+h) 

Ci,i(A) = Ci,iA- a,N+38 n (A i + A)(A l -A), (4.25) 

i=i 

2lN-(a+b) 

Ci, 2 (A) = Cl)2 ^+6+N A -aN+2a JJ (Aj + A)(6jAj - Xq), (4.26) 

1=1 

where the e$ = ±1. The property (1C.12K the relations (I4.18t and (14.171 imply that Vi, e; = 1. Let us now introduce the poly- 
nomials Ci,i(A), Ci, 2 / + i(A) and Ci j2 (A) defined by simplifying the common factors in Ci,i(A)/Ci i i, Ci j2 j-|_i(A)/Ci )2 ; + i 
and Ci, 2 (A)/Ci j2 , respectively. Then, Ci.i(A) has the form: 

Ni.i 

Ci,i(A) = J\(Xh - A), with Ni 4 < 2ZN - (a + b). (4.27) 



h=l 



°The cii.j and bij are non-negative integers and 7^ for any h G {1, 4ZN — (aij + bij)} 
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Here and in the following, we will write Zf(X) to denote the set of the zeros of the function /(A), then it holds: 

^Ci,i(A) c — ) A 2 iN-(a+6)}. °^Z Cii(A) , S P) A <t tl (A)> (4.28) 

for any S p .a = (Ao 7 9A0, g 2 'A ), a p-string of center A G C. Now, since by definition n ,y. = Z^ n 

Z^ 2(+1 (a) = 0, equation ( I4.171 i implies: 

Ci, 2 /+i(A) = ^-^nCAg- 1 ), C 1)2 (A) = q- N ^C 1A (Xq). (4.29) 

Then, equation ( 14.201 can be written as a Baxter equation for the polynomial Ci t i (A): 

i(A)Ci,i(A) - a(A)Ci,i(Ag- 1 ) + d(A)C 1 , 1 (Ag), (4.30) 

with coefficients a(A) = q Nl l tp&(\) and d(A) = q~ Nl - 1 (p~ 1 d(X) and ip ^C^i/Ci^ =Ci. 2 ;+i/Ci ! i. Note that the consis- 
tency of the above equation implies that ip is a p-root of the unit and then, by (14.23K we have tp = q@ +Nl \ 

Let us define now the polynomial 

2/N-(o t +6t) 

Qt(X) = A at Cij(A) = A a * [J (X h -X), (4.31) 

h=l 

where: 

Ni,i = 2iN-(ot+6t), a t G {0, 2Z}, b t = b + mp, m G Z-°; (4.32) 

then g~ at = g Nl l (p, Qt{X) is the desired solution of the Baxter equation (14. 1U which belongs to C[A] 2 ;n- The character- 
izations ( 14.121 ) are then trivial consequences of the asymptotics of f (A) € plus those of the coefficients of the Baxter 
equation ( 14. 1 lb satisfied by Qt(X). □ 

Lemma 7. The polynomial solution Qt{X) of the Baxter equation \4.11\ associated to t(A) G S T2 ;i a (-real polynomial, 
i.e. it satisfies the following complex-conjugation condition: 

(Qt(A))* = Qt(eX*) VA e C, (4.33) 

where e = ±1 is the discrete parameter defined in ( 12.251 ). Then, the Baxter Q-operator, Q(X) defined in ( 14.71 ), is (.-self- 
adjoint and it is a polynomial of maximal degree 2/N in the spectral parameter A. 

Proof. Let us remark that the self-adjointness of the transfer matrix r 2 (A) and its parity properties imply: 

(t(A))* = t(X*) and t(-X) = (-lft(X) VA 6 C, (4.34) 
while the coefficients of the Baxter equations satisfy the complex-conjugation conditions: 

(a(A))* = e N d(eA*) VA G C. (4.35) 
Then, under complex conjugation, the Baxter equation reads: 

i(eA*) (Q t (A))* = d(eA*)(Q t (A<r 1 ))* + a(eA*)(Q t (A (? ))* VA G C, (4.36) 

and so the injectivity of the map t(A) G E T2 — > Qt(X) G C[A] 2 ;n implies the condition (14.331 1. □ 

Remark 4. For any t(X) G £ r2 the Theorem|4]defines a link between the corresponding polynomial solution of the Baxter 
equation and a determinant of a (p — 1) x (jp — 1) tridiagonal matrix. Let us note then in literature there exist also other 
determinant characterizations of solutions of the Baxter equation. In the quantum periodic Toda chain, linear combinations 
of determinants of semi-infinite tridiagonal matrices allow to express these solutions [61, 62, 63 1. A careful analysis of the 
r 2 -model in the limit /3 2 — > p 1 with p 2 irrational (i.e. p', p — > +oo) is of clear interest as the dimension of the representation 



2 'The proof of this statement coincides step by step with the proof given in Theorem 2 of iTfl. 
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of the Weyl algebra W n and the size of the tridiagonal matrix, related to the Baxter equation solution, diverge. Subsequently, 
one could ponder whether the characterization which holds in the case of the quantum periodic Toda chain also applies to the 
T2 -model for irrational /3 2 . This could lead to a reformulatior@ of the t-i -spectrum in terms of solutions of nonlinear integral 
equations (NLIE) related to the results in ll66l . Interestingly, another DDV-typ^| NLIE reformulation can be introduced 
and it describe the complete T2 -spectrum thanks to the completeness of the Bethe ansatz. 

4.3 Completeness of Bethe ansatz type equations 

The previous analysis leads to the complete characterization of the transfer matrix spectrum (eigenvalue and eigenstate) in 
terms of real (in the case e = 1) or imaginary (e = —1) polynomial solutions of the Baxter equation (14. lit . 



Proposition 2. To any T2-self-adjoint representation which further satisfies the conditions \4. 3 i is associated a system of 
Bethe ansatz type equations and there exists a one to one map between the set of all the p-strings free and e-self-adjoint 
solutions of it and the set E T2 . 

Proof. Let us construct the isomorphism. Theorem|4]and the previous lemma imply that to any t(A) G £j: 2 it is associated 
one and only one p-strings free and e-self-adjoin@ solution (Ai, . . . , \2iN-(a t +b t )) of the system of Bethe ansatz equations: 

, v 2JN-(at+6t) , , 

SH 2 "' n Bf^tr * e {i,...»N- («.+*», <« 7 > 

with a t — ±fc mod p, b t = ±fc mod p. 

Conversely, given a p-strings free and e-self-adjoint solution of ( 14.37b with a = ±b modp, then it defines uniquely a e-real 
polynomial Q(X) by the equation ( 14.10b . Now, by using Q(A) we can construct the function: 

t(X) = (a(A)g(< Z - 1 A)+d(A)Q(gA))/Q(A), (4.38) 

which thanks to the Bethe equations ( 14.371 ) is nonsingular for A = A c , Vc S {1, . . . , 2ZN — (a + b)} and is a e-real Laurent 
polynomial of degree N in A. Moreover, we can also construct two states | t± a ) by inserting Q(A) in equation ( 13.471 ). So, 
t(X) and ^t(v) = {Vli 1 1 ) satisfy the system ( 13.15b and ( 13.16b . for k = ±a. Then the fixed solution of the Bethe 

equation allows us to reconstruct uniquely the T2-eigenvalue t(A) and the simultaneous eigenstates | t± a ) of T2(A) and 6, 
with 0-eigenvalue q ±a . □ 

Note that the previous results can be rephrased as completeness of the Bethe ansatz type equations generated by SOV, as 
from this type of Bethe solutions we can reconstruct the complete set of T2 -eigenstates. 



5. Spectrum characterization of the inhomogeneous chiral Potts model 

In O Bazhanov and Stroganov have proven the following remarkable connections between the integrable chiral Potts model 
and the -model: 

i) The fundamental R-matrix intertwining the T2-Lax operator in the quantum space is given by the product of four chiral 
Potts Boltzmann weights. 

22 NLIE has been used first to reformulate the spectrum of integrable quantum models in 1641 1651 . Similar NLIE were presented in 1 66 1 
and also in (671Pl[55lp). 

23 This type of NLIE was introduced and analyzed in |70 |-[77 1 for the fermionic lattice regularizations of the sine-Gordon model, see 
also 1781 - 18 II for a related model. 

24 It means that the tuple (Ai, . . . , XziN-iat+bt)) satisfies the following property under complex conjugation: 

(Al, . . . , A 2 JN-(oi+6 t )) = ( e ^ir(l)> • • • ! e -^7r(2!N-(a t +i>t))) 

where n is a permutation of the indexes {1, 2/N — (at + bt)}- 
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ii) The transfer matrix of the chiral Potts model is a Baxter Q-operator for the T2-model. 

Both these statements are true when the t 2 -model is restricted to parametrizafion by points on the algebraic curves^ C^. In 
this section, we use the property ii) to characterize the eigenstates of the inhomogeneous chiral Potts transfer matrix by the 
SOV construction. 



5.1 Transfer matrix of chiral Potts model as Baxter operator of T2-model on the chP curves 

Let us consider the following tuple p= (a p , b p , c p , dp) £ C 4 and let us introduce the following notations: 

x p = a p /d p , y p = b p /c p , s p = d p /c p , t p = x p y p , k 2 + (k ) 2 = 1, (5.1) 
then the algebraic curve Ck of modulus k is by definition the locus of the points p which satisfy the equations: 

xP + yP^k(l + X PyP), kx% = l-k's-P, kyP = l-k' S P. (5.2) 
Similarly to [3|, we can introduce the parametrization: 



Xa n = -t p 1/2 b qn b rn , hn/rq 1 / 2 = (x p /y p ) 1/2 a qti d ln /q 2 , 

/3 n /A = -4 /2 rfq„rf r „, q 1/2 a, n /r = -~(x p /y p ) 1/2 c qn b [n , 

S n X = q- 2 t p 1/2 a qn a In , q 1/2 d n r = -(y p / x^^d^a^, 

7„/A ee tj /2 c q?i c r?l , q~ 1/2 <£ n r = (y p /x p ) 1/2 b qn C ln , 

of the T2-Lax operator in terms of the points p, q„ and r„ of the curve (15.21 i of modulus k. Note that we can identify: 



(5.3) 



— 1/2 

so that the off-diagonal elements of the T2-Lax operator do not depend on the point p on the curve and the parameter t p 



r ee {y p /x p ) 1 ' 2 , (5.4) 

. open 

is then proportional to our spectral parameter: 

tp- 1 / 2 ee A/c . (5.5) 
Then, we can introduce the operator T^ hp , characterized by the kerneQ 

N 

Tf (z,z') = < z |TfV> = J] W qnV (z n /z' n )W tnP (z n /z' n+1 ), (5.6) 

n=i 

in the left and right u„-eigenbasis (12.5b . where the functions W and W are defined in ( IB. 11 1. T^ hp coincides with the transfer 
matrix of the inhomogeneous chiral Potts model 0. As anticipated, T^ hp is a Baxter Q-operator w.r.t. the T2 -transfer matrix: 



where the coefficients read: 



T 2 (A)Tf = a BS (\)Tf? q + d B5 (\)Tf*, (5.7) 



N 1 | <? 1/2 Ab„ 

i ^ Priori 1 H Tjo 

n— 1 q ' r<c n 

N , i _ Ad„r 

*s(A) = (-ifYlM^jl + ^X) (5-9) 

n=l 1 n n a „ 

W qnP (z(l)) W rnP (z(l)) 

^-w qnP (z(0))w ln M0)Y ( 3 



and: 



25 See the next subsection for precise definitions. 



26 For a direct comparison see formula (4.12) of 01221 with the following identifications: 

Zj =q 2 ^, z' 3 =q 2 ^ VjG{l,...,N}. 
Note that T c ^ p is well defined, the W-functions dB.lt being cyclic functions of their arguments. 
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Note that to write these coefficients we have inverted the formulae (15.3b : 

x <\ n 3/2 Ab ra x In _ 1/2 Ad„r csm 

2/p Pn?* Pn 

2/q„ _ „-i/ 2 ?'Aa» y r „__ i /2 Aa„ « 

— U i — H i 6 q„*r„ — • W- iz 7 

It is worth pointing out that while the Baxter equation (15.7b holds in the general inhomogeneous representations, the com- 
mutativity properties: 

b2(A),Tf]=0, [Tf\Tf]=0 VA^eC (5.13) 
only hold under the further restrictions: 

q„=r„ Vne{l,...,N}. (5.14) 

Moreover, we also have: 

as it simply follows by computing the matrix elements on the left and right u„-eigenbasis: 

N 



e,Tf]=0, (5.15) 



(zieTfV) = ( Zl /q, z N /q\Tf p \z') = J] W qnP (z n /qz' n )W TnV (z n /qz' n+1 ) 

n=l 

= (z\Tf\qz'- L ,...,qz' N ) = (z\TfQ\z'). (5.16) 



5.2 General representations: spectrum simplicity and eigenstates characterization 

Let us denote with 1Z C ^ P the subvariety of the r 2 -representations of real dimension 4N + 4 parametrized by: 

a n = -&q n /c , ton = -dn/q = -a^d^/ q 3/2 , (5.17) 
Pn = -cod, , c„ = -Si n q = & qn c qii q 1/2 , (5.18) 



am _ 

c eC, q„eC k , keC. (5.19) 

Now, for the representations 7\LjJ p , we can prove: 

Proposition 3. All right and left eigenstates of the transfer matrix T A are eigenstates ofT2{X) with corresponding eigen- 
values related by the functional Baxter equation: 

_ mk^w ( , 20) 

where Cf£ p denotes the T C £ F -eigenvalue. Moreover, for almost all the representations in TZ L ^ P , Cf£ p is non-degenerate and 
the simultaneous eigenstate of(T c £ p , r 2 (A)) is the one associated to the eigenvalue t(A) as constructed in Theorem\3\ 

Proof. The first statement concerning right eigenstates is a simple and well known consequence of the Baxter equation 
(15.7b from which we also have that the corresponding r 2 -eigenvalue t(X) is given by (15.20b . To prove the statement for 
left eigenstates we also have to use the first commutation relation in ( 15.13b . Now, the fact that the r 2 -spectrum is simple 
for almost all the representations in 7?.* p implies the remaining statements of the proposition, once we recall that the r 2 - 
spectrum characterization of Theorem [3]holds for completely general representations and so in particular for 7^j hp . □ 

Remark 5. There exists a proper subvariety of 1Z C ^ P for which t 2 (A) has double degeneracy eigenvalues. In these repre- 
sentations we have shown that the degeneracy is resolved by the charge O then, thanks to the commutativity ( 15.151 ). it is the 
simultaneous spectrum of and T A hp that is simple here. 

27 Note that the dimension of TZ^ r is defined by counting two independent complex parameters for each points q ?l £ Ck- For example 
o qri and d qn can be taken as independent complex parameters while the other two b qji and c q?i are fixed by the equations of the curve Ck- 



22 



5.3 Self-adjoint representations: complete spectrum characterization 

The Proposition [3] characterizes completely the T^ hp -eigenspace showing that it is contained in the T2-eigenspace. The 
reverse inclusion of eigenspaces however is not proven for general IZ^. To clarify this point it is worth to remark that while 
T^ hp and T2(A) are both one-parameter families of commuting operators, their mutual commutativity is proven only when 
the spectral parameters coincide, see (15.13) . This means that we cannot argue, by using this commutativity onhl^l that the 
T2-eigenstates are eigenstates of T^ hp ; in fact, there may be some residue degeneracy in the T2-spectrum for fixed value of the 
spectral parameter. Note that the simplicity property of r-i -spectrum only says that, taken any two different T2 -eigenstates, 
the corresponding two T2-eigenvalue functions are not identical. This is a less strong statement that: for any T2-eigenstate 
t ) there exists at least one value At of the spectral parameter such that the T2 -eigenvalue in X t is non-degenerate. If this 
last statement can be shown, the commutativity ( 15.13b implies the coincidence of the and T* p -eigenspaces. 

Of course, the diagonalizability of T^ hp defines a sufficient criterion to overcome deeper analysis of the T2 -spectrum degen- 
eracy. Indeed, in this case the Proposition|3]directly implies that the T|J p -eigenbasis is simultaneously a T2-eigenbasis. Then 
it is relevant to know the precise characterization of the representations of 1Z C ^ P for which the chiral Potts transfer matrix 
T^ hp is proven to be diagonalizable. Here, we prove that this holds in particular for the representations on the non-trivial 
curves for which the -transfer matrix is self-adjoint. Let us denote with 7^ hp - s ~ ad j = TZ^ V D IZ^^ such subset of IZ^, 
then the following characterization holds: 



Lemma 8. The subvariety j^ p - s ad j [ s characterized by A5.17i - i5.19i and the following points q n on the chP curves: 

In = ( a 9ni e 1 e o,nag n ,e ,ndg n ,d qn ) e C k , e ,„ = ±1, k* = efc, (5.21) 
where we have fixed Cq = 1, just to simplify the parametrization. 

Proof. This subvariety is of real dimension 3N + 1 and it is parametrized by: 

( a n) 2 

Pn = IL ^, d„ = ~ea* , c„ = -a„g, b„ = ea*/g, (5.22) 

with: 

a n = -(qa*J 2 , a„ = -eq 1/2 a* n d* n , (5.23) 

where the point q„ is characterized by (15. 2U . Indeed, it is a trivial check to verify that the above parametrization satisfies 
both the condition ( 12.251 ) and those defining ^ hp . The e-reality condition on the modulus k of the chP curves is then derived 
from the equation of the curv 

k = , -4P , , (5.24) 

where we have denoted x q = ip Xq \x qn \ with (p Xq the phase of x q?i . □ 

Lemma 9. Let q be a point on the algebraic curve characterized by the condition ( 15.271 ), then the W-functions ( IB. 71 ) satisfy 
the following property under complex conjugation: 

W,Mn))y_W qM P-n)) yz£ ^ (525) 



W qp (z(0))J W qp (z(0)) 
where we have the following restriction on the point /?£ Cfc, with k* — efc: 

(5.26) 



28 Note that the reverse statement is proven thanks to the Baxter equation d5.7> 



9 It is worth noticing that for e = 1 the equation J5.24b implies the following restriction on the modulus of the curve fc e] — 1, 1[. 
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Proof. It is important to point out that the restrictions on the point p are compatible with the condition p£ Cfc. Now, by the 
definition ( IB. 11 ) and the condition ( 15.2U satisfied by the point q and (|5.26l l, we have: 

W *M™)) \ * = (ss) -n tt y P -i 2k y* (5 27) 

where we have used that s* = s^ 1 and s* = s p . Once we use the cyclicity condition satisfied by the W qp (z), the r.h.s. of 
the above equation coincides with the r.h.s. of ( 15.25b . □ 

Proposition 4. For any representation in lZ c ^ P ' s ad \ the inhomogeneous chiral Potts transfer matrix T°£ p is a one-parameter 
family of normal operators w.r.t. the points p defined in ( I5.26D . In particular, it satisfies the following Hermitian conjugation 
property: 



A %y*(o))w?>(Q)) 



-rchP\ \ _ „ xolP 

, A l(P,?„) J ~ 9 <P,4 n ) A 



w/zere T^ p is the second chiral Potts transfer matrix defined by: 

N 

(z,z') = (z|ff |z') = 1] ^(^i/O^wWO- (5-29) 



Ta 
X 



n— 1 



Proof. The proof of (15.281 1 is a simple consequence of the previous lemma, then the known commutativity 1 3 1 : 

w»4 n ; (p'.q n ) 

implies the statement of normality thanks to ( 15.26b . □ 
Let us denote with S T d,p the set of all the eigenvalue functions of T^ hp , then: 

Theorem 5. For any representation in 1Z C ^ P,S ad ' the formula ( 15.201 ) defines a one-to-one map between the eigenvalue sets 
S T t*p and £ T2 . Moreover, Theorem\3\allows to construct the full (T x hp , T2(A), Q)-eigenbasis associating to any t(X) 6 S T2 
the corresponding eigenstate. 

Proof. The diagonalizability of T* p implies, thanks to Proposition [3] that T^ shares with T2(A) and 9 a complete set of 
eigenstates from which the statement about the eigenvalue functions follows. Then, the map ( 15.201 ) being an isomorphism, 
we can label all the eigenstates of T* p in terms of the eigenvalues of r 2 (A) and 9 and, by using Theorem[3] we can construct 
the entire simultaneous eigenbasis. □ 

5.4 Completeness of Bethe ansatz type equations 

Let us denote with 7?.^ hp ' wd ^ the subset of 7£j* p > s " ad j which contains the representations which satisfy the constrains ( 14.31 ). 
Then the following characterization holds: 

Lemma 10. In the case o/N even, the subvariety TZ C ^ P "' 1 is characterized by f l5.77D - 0.79l ) and the points q n on the curves 
f tOTT l with moreover: 



q2n+i = ±elil i <p2- 2S ^<? s *-°\a q2n J, d q2n+i = ±e^ n+i e^t 4 K 2 „ +i |, * = 0, 1, (5.31) 



where ip aq is the phase of a qi and 

\\ £2,h, ei,n = ±L £2,n = ±1- (5.32) 



n-l 

e n = 

h=l 



In the case of N odd, lZ c ^ ps ad ' is obtained by the further requirements: 



V l=±q 2 , ei,„ = e Vn G { 1 , . . . , N } . (5.33) 
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Proof. Let us impose the second and third conditions of (14.3b . these together with (15.22b imply: 

<y=^ <±i = ^ ; (5 34) 

a„ / a„+i a* 

from which the first condition in (14.3b is automatically satisfied for N even while it imposes a n G M for N odd. The 
conditions ( 15.34b imply for the points q„ the constrains: 

rf q„ „2 a q„ r 2 a q. + i _ _ a q. 
d q„ °q„ fl q„+i 1 fl q n 

with ei.„ = ±1. Then, the parametrization for the points q„ is given by ( 15.21b plus ( 15.31b and the equations of the curves 
for these points q„ read: 

In the case N even we do not have to add further constraints, and it is then clear that the subvariety 7£jjj lP ' s J can be for 
example parameterized by the N + 2 parameters (k, \d qi \, \a qi \, \a q \, \a qfi |). In the case N odd we have still to impose 
a n G R, which reads: 



In this case, only the choices: 



lead to non-trivial curves, k ^ 0: 



< = ±q 2 - (5-37) 



ei,„ = e VnG {1,...,N}, (5.38) 

±26„e 1 /2 



eKJf + 1/K|f 



(5.39) 



and so it is clear that the subvariety 7^ hP ' s ~ adj can be for example parameterized by the N + 1 parameters (k, \a qi |, |a q J). 

□ 



The interest toward these representations is due to the complete characterization of the T2-eigenbasis by the e-self-adjoint 
solutions of the associated Bethe equations ( 14.37b . 

Theorem 6. For any representation in 1Z C ^ P ' S ad \ there is a one to one map between T, T d,p and the set of all the e-self-adjoint 
solutions of the system of Bethe ansatz type equations A4.37i . Moreover, the two simultaneous eigenstates ofT c £ p and O 
corresponding to the given (f* p G E-pw are then constructed by inserting in equation A3. 47]) the e-real polynomial Q(X) 
defined by \4.10b in terms of this Bethe equation solution. 



Proof. The isomorphism stated in the above theorem is just the composition of two isomorphisms. First, since by definition 
7£.chP,s-adj c 7\L^ hP ' s " adj , the isomorphism defined in Theorem |5]holds and (15. 20b associates to any eigenvalue q^ hp £ S T ,±p one 
and only one <(A) £ S T2 , and viceversa. Then, we can use the Theorem|4]to uniquely associate the e-real polynomial Q(X) 
defined by ( 14.10b whose zeros define a e-self-adjoint solution of the Bethe system of equations ( 14.37b . The uniqueness of 
the polynomial Q(X) is proven as in Theorem[3] indeed for any representation in 7^ hP ' s '~ ad \ the coefficients of the Baxter 
equation ( 14.11b have different average values: 

m^m = ( _ irT r (^KiA) 2 



□ 
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6. Outlook 

One of our motivations to the present work is to define the required SOV setup to generalize the approach used in [82] to the 
T2-model associated to the most general cyclic representations of the 6-vertex Yang-Baxter algebra, such a generalization 
will be presented in [ 121 1. There, the obtained reconstruction of local operators in terms of quantum separate variables and 
the derived determinant formula for the scalar products of states are used to compute the form factors of local operators 
on the transfer matrix eigenstates and to express them as sums of determinants given by simple modifications of the scalar 
product ones. 

It is worth mentioning that this approach is not restricted to cyclic representations of 6-vertex Yang-Baxter algebras, as it is 
addressed to the large class of integrable quantum models whose spectrum (eigenvalues & eigenstates) can be determined by 
implementing Sklyanin's quantum separation of variables. Such a statement has already been verified for some key quantum 
integrable models associated to highest weight representations of the Yang-Baxter algebras and their generalization; these 
results will appear shortly in a series of papers. 
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A. Constructive proof of the existence of cyclic SOV representations 

In the following subsections we will show by recursive construction the TheoremQ] 
A.l Recursive construction of B-eigenstates 

Similarly to [1|, we will construct the eigenstates (77 1 of B(A) = B N (A) recursively by induction on N. In particular, 
assuming we have the B M -eigenbasis for any M < N, the eigenstates { 77 1 of B N (A) can be written in the following form 

^1 = EE^nMxJxJUJ^xJ, (A.i) 

in terms of ( x z I an d ( Xx l> respectively, the B M -eigenstates and the B N _ M -eigenstates with eigenvalues defined as in ( 13. U 
by the tuples Xz = (Xzo)o=i,...,m and Xx = (Xxo)o=i,...,n-m- Computing the kernel K N (r)\ x z \x x ) will give the recursive 
construction of the states ( 77 1 . 

SOV- representation for N = 1: let us construct first the SOV representation for N = 1. We introduce the states: 

/ (0) h t I _ 1 llr=lW ' ^1 + q v 7 ; lt>i) 

(>7lV| = 2^ (a p + foP) fc /p ( ' VAl6 ^ * (A ' 2) 

where (1, k\ are the states of the Vi-eigenbasis. ' q hl | for hi S {1, ■■■,p} then defines a Bi-eigenbasis: 

(4° ) <l hx \Zi=vi 0) q hl {Vi 0) q hl \ with (r# )) ) P = g" /a (af + (A.3) 
Subsequently, Ai(A) and Di(A) have the following representation in the Bi-eigenbasis: 

(»?1 V 1 |A 1 (A) = Aax^V^VftA-^V^ 1 !, (^V 1 |Di(A) = 11 X- 1 ( V f ) q h ^ 1 \-d 1 X( v f'>q h ^ 1 \. (A.4) 

\A.l\ a Reconstruction of the kernel K N (r] | X 2 |Xi ) w.r.t. \z o=1,...,m-i and Xi o=i,...,n-m-1- 
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The decomposition 

B N (A) = A 2 M (A)B 1 n _ m (A) + B z M (A)D X N _ M (A) (A.5) 
implies that the matrix elements of the kernel K^(rj \Xz\Xi) satisfy the equations 

\ n Z M 

(A)B 

1 N — M 

(A) + B 2M (A)D 

1 N — M 



N-l 

V\Xz,Xi): 



(A.6) 



a=l 



where stays for the transpose of an operator 0. Taking 

XLc.q hl $ Z dGtqMi n _ m(A ), Xzbq h ' 2 i Z dotqMi m(a) and xi a q hl j= Xzbq h \ (A.7) 

where hi e {1, ...,p}, a S {1, N — M} and b G {1, M}, the previous equations can be written as recursion relations 

for the kernel w.r.t. Xia an d Xzb by simply fixing A = Xia=i n-m-i an d A = X26=i,...,m-i- Indeed, for A = Xia the 

( IA.6I ) leads to 

M-l 



K Av\xz\q Sa x x ) d i sov) (q 1 x la ) Xzu Y[(XxJxzb-Xzb/Xi a ) 

b=l 
N-l 

= k n{v\Xz;x3.) % n(x la /Vb-Vb/x la ), 



(A.8) 



(A.9) 



b=l 

while for A = Xza it holds 

N-M-l 

K Jv\q Sa Xz\X L ) a. ( z SOV) (q +1 X za ) XxN-M J| (^a/Xxb -Xxb/Xza) 

fc=l 

N-l 
b=l 

Here d[ SOl/ ' 1 (x la ) and &z S ° V \xza) ^ the known coefficients of the SOV representations in the subchains l and 2 while 
q ±s - is defined in dTTTb . 

IA.7I £> Reconstruction of the kernel K N (r] \ XzlXi ) w.r.t. Xz m fln ^ X± n-m- 

In the T2 -model the form of the asymptotics of the Yang-Baxter generators w.r.t. A leads to some complication in the 
computation of the recursions satisfied by the kernel w.r.t. the variables x z u m ^ Xi N - M - H ere > we show that just doing a 
discrete Fourier transform w.r.t. these variables we get a simpler characterization. Let us introduce the state: 



(vPq h «,r, i °±iq hN - 1 ,-,ri{ 0) q hl \ = E ^ (vi°^ kn ,Aq h "~\ ---A^^l (A. 10) 

fc N =i 



then it is an eigenstates of the 8-charge: 



(v { ° ) q h ",{vi}\®N = q h "(riN ) q h ™,{Vi}\ (A.ii) 

and moreover it holds: 

uii V 7 ?" A / 

Note that by using the recursion formula (1A.U and the representation ( lA.lOb we get: 

(^VM^}! =J2J2 R MXz\Xx){Xz, m,Xz, a I ® (Xx, n-m,Xi, al, (A.13) 

with 

KMxz, M q h ^ M ,{xz, allxx, N - M g^ N - M ,{xx, a}) = 
1 p 

V g-(kx. N-Mhx, N-M+*». Mfcx, M)^^^ M q k ^^,{x z . a}|Xx, N -M 9^ ' ^ " , {Xx, a})- (A. 14) 

^— ' _ 

&2., Mi^i, N-M — 1 
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To obtain the dependence of K N w.r.t. \z, u an d Xx, n-m> let us introduce the asymptotic operators Bj, whose action 
reads: 

N-l 

(r,\B^ = ( T 1)^\ Urn A ± ( N - 1 )(r / |B N (A)= % n^ ±1 (^L (A.15) 

log A— >=poo 

z=l 



while the decomposition ( IA.5t implies: 

N-M-l M-l 

Bl T) - a^efxx, n — m I] *M + m n *M- (A. 16) 

i=l i=l 

Acting on the generic state (t/n°\ {?7i}| with we get a system of two equations: 

KM XzlXx ) = *e^ N (»7 | Xzk" 5N " M Xx ) + Y^nU I <T* M XzlXx ), e - ±1- (A.17) 
Here the coefficients read: 

a €,Z<I Al, N-Mllj = l Xl, i "e,l ( / Xz, Mllj=l Xz, i .. 

X ^ = t-tN-1 - e . = n N-l - t ■ ( A - 18 ) 

Then solving this system we get our result: 

K^vlq'^xJXr) x_-x+ K^rilxJq-^^Xx) Y+ - Y- 



Finally, ( |A.20t leads to the following recursion relation which fix the dependence of the kernel w.r.t. %: 



km»v i xz; xx ) _ xr^xrd m=T x 2 fc njTi x. 



*U v I <z 5m x 2 ; <z 5n - m Xx ) 77^ +) n^i 1 % 



A.2 Gauge-invariant SOV dates: Z r , Z { ^\ Z^\ A{Z r ), V(Z r ) 

The recursion relations (lA.8t - (lA.9) and (1A.19) and the requirement of cyclicity give a system of N algebraic equations in 
the N unknowns Z a =r\ v a : 

M-l N-l 

^xN-m(xL) Xzm II (Xla/Xl b ~ Xljxla) = V> J] ^IJ< ~ Vb/xL) , (A.25) 

6=1 6=1 

N-M N-l 
^zm(xL) Xxk-m II (X5a/X? 6 - xVxSJ = T& R ixlaK ~ </xL) , (A-26) 



(A.19) 



^nC^IXzIXx) x_y + -x+y_ # N OIXzlXx) x_y + - x+y_ 

\A.l\ c Reconstruction of the kernel K N ( rj \ x 2 ; Xx ) w.r.t. r\. 

The formula 

D N (A) = D 2 M (A)D 1 n— m(A) + C 2 M (A)Bx N _ M (A) (A.20) 

implies the identity 

<77|D N (r7i) = -det q M 2 M (qr/ l )(r/\B^ 1 M (qr/i)B 1 N _ M (l7i), (A.21) 
writing C 2 M (A) by the quantum determinant det q M 2 M (A) in the subchain 2 and using 

(r?|Bi N _ M (?7i) = -(^(Dj. N - M (f?i)A~ ^(g?7i)B 2 M (<?»?*). (A.22) 

Then, the following recursion relations is obtained for kernel if N ( ?7 1 x 2 ! Xx ) w - r -t- »/t=i,...,N-i: 

K N (q s "i] | X 2 ;Xi ) _ det q M 2M (q^) 

XlN — M n a =l {m/Xia-Xiahi) 



(A.23) 



(A.24) 



6=1 6=1 
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and 

N-l N-M-l M-l 

zn n z t x = (±i) n ^«i ±z, n m n z i« + (i 1 )"" 1 ^, ±^ m n z ti v&n 

n—l a—1 a—1 

where we have used (13.6b in the subchains 1 and 2 to express Z^\ Z~^\ and Z^q Ila^i 11 ■^a- Th e above system 

of equations completely determines the Z a in terms of Xza, X?a, as we can rewrite it in terms of the following Laurent 
polynomial equation: 

[N] 

Ai(A)£ N _ M (A) + S m (A)P n _m(A) = Z n [] (A/^ a - Z a /A), (A.28) 

a=l 

as a consequence of the simplicity of the spectrum of B(A) in both the subchains 1 and 2. Note that the l.h.s. of the above 
equation is formed out of the known average values of the monodromy matrix elements in the subchains 1 and 2. So the Z a 
for a 6 {1, N — 1} are fixed determining the zeros of the known Laurent polynomial at the l.h.s. of ( |A.28t . The cyclicity 
allows to fix the remaining gauge-invariant SOV dates: 

V H {Zi) = -dct q A^ZO g "~^y , AuiZi) = -det q A* N _ M (Zi) J^ftL (A.29) 
for any i G {1, N — l} where: 

v 

det q 7Wx(A) = IJdet q M x (<fA), X = M,N-M,N. (A.30) 

a=l 

where it holds: 

det q MxW = Ax{A)V x (k)-Bx(A)Cx(A), X = M,N-M. (A.31) 
A.3 B-spectrum completeness and simplicity 

Let us prove that the previously constructed set of B-eigenstates ( r\ | is complete by showing that there are p N distinct 
B-eigenvalues. 

Proposition 5. The SOV dates Z r with r G {1, N — 1} are all distinct for almost all the values of the parameters 
oP n , ft?, a£, cfj, of the T2-model. 



Proof. From the functional dependenc^j of Z x , . . . , Z N1 w.r.t. these parameters we just have to prove that the Jacobian: 

J(a p , (3 p , a p , b p , c p , d p ) = det ( ) ^ (A.32) 

does not vanish for some special values to derive J ^ for almost all the values of the parameters. In fact, for J ^ 0, the 
map Z = Z(a p , . . . , o^-i) * s invertible from which the claim of the proposition follows. 

To show that ( 1A.321 >. let us consider the representations which satisfy the following constrains: 

b£+a£ = 0, < + d£=0 Vn € {1, N - 1}, (A.33) 
while we leave the representation in the site N unrestricted. Then the averages (13.121 simplify to: 

A <~^ A ° ) Vne {!,..., N-l} (A.34) 

Then by ( 13.111 ) we can compute the l.h.s. of ( IA.28b and we obtain: 



<l 



P/2 



« + to n ^<^r - S» = Zn n<i - f >■ 

and then J / trivially follows. We can then define r/ a °^ as a p-root of Z a . They satisfy ^1°^ 7^ ^ or an y 

a ^ b e {1,...,N - 1}. □ 
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Finally, keeping in any subchain the representation of the B-spectrum simple, the Theorem[T]follows by induction. 



B. Baxter Q-operator construction for general r 2 -model 



In ifTTl . Baxter has extended the construction of the Q-operator by gauge transformations of [3 1 to the T2-model for general 
cyclic representations not restricted to those parameterized by points on chP algebraic curves. The main tool there was 
a generalization of the discrete dilogarithm functions. In particular, he has remarked that asking the cyclicity only for 
the products of couples of these functions does not prevent the Q-operator from being well defined for general cyclic 
representations of the model. Here, we reproduce and adapt Baxter construction in our notations and we point out some 
interesting connection with the SOV construction. 



B.l Dilogarithm functions on the algebraic curves Ck 

B.l.l Dilogarithm functions 

The dilogarithm functions have been previously introduced and analyzed in [45 , 55 1. Here we use the notation: 

where z(n) — q~ 2n , Vn G {0, 21}. They are solutions of the following recursion relations: 

W qv (zq) _ z s v x p ^-^q*- 1 W qp (zq) _ qz' 1 y p 1 ~ 
Wqpizq- 1 ) s q y p l-^-q-iz' W^zq- 1 ) s p s q x p 1 - ^q~ 1 Z - 1 ' 



If the points p and q belong to the curves Ck, they are well defined functions of z G § p which satisfy the cyclicity condition: 
fi.i.2 Recursion relations 

It is worth defining the following three discrete automorphisms of the Ck -curve: 

i) 5 : p = (a p ,fep,c p ,dp) G C k -> H(p) = (ga p , <?6 P , c p , d p ) G C fe , (B.4) 
m) A : p = (a p ,6 p ,c p ,d p ) G C fc -> A(p) = (gfe p , a p , gd p , c p ) G C k , (B.5) 
m) T : p = (a p ,6 p ,c p ,dp) G Cfc T(p) = (-g 2 dp,c p ,6 p ,-a p ) G Cfc. (B.6) 

If the points p and q belong to the curves Ck, the ly-functions satisfy the following recursion w.r.t. the action of 5 on the 
point p: 

w qp {gz) = 1/a jM£(0) l - x ji~ 1 w qp {q Z ) 1/2 iy qp (z(p) i - f^' 1 

WqBCp)^) Z ^(0)) 1 - |<rV V^jW ^ ^qp(^(0)) ' 

and 



Wgp(g*) = -1/2 gWQ) 1 - t q ^ _W qp {qz)_ = _ 1/2 W qp (z(l)) 

WqH(p)(z) Z VK qp (z(0)) 1 - |g-l ' W- qS -i(p)(«) Z W'qpWQ))!-^- 1 ' 

B.2 Parametrization of r 2 -model by points in C 3 

Let us implement the following gauge transformation on the T2-Lax operator: 

L ™ A V l/r n+1 z' n+1 1 J L " A i -l/r„< 1 J' 
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where r n £ C and 

z' n eS p = {q 2h , h = l,..,p}. (B.10) 
We associate to the tuple (a v , b v , c v , d v ) € C 4 the following point: 

v = (x v = a v /d v , y w = b v /c v , s v = d w /c w ) £ C 3 , (B.ll) 

and we define: 

t v = x v y v . (B.12) 
Now, let us consider the point p = (xp, z/p, Sp) £ C 3 and from it the points: 

P„ = (Xp/Vn, VpVn, Sp(T„) £ C 3 , O n € C, V?l £ {1, N}, (B.13) 

then it holds: 

tp n =t ¥ . (B.14) 

We can introduce the following parametrization: 

Xa n = ~t ¥ 1/2 b qn b Vn , hn/rnq 1 / 2 = (Zp n /2/p-J 1/2 a q „£k„/<? 2 , 

/3 n /\ = -tl /2 d qn d rnl q 1/2 a, n /r n = -{xp n /ypj 1/2 c qn b lnl (B 15) 

S n X = q~ 2 tp 1/2 a qn a In , q 1 / 2 d n r n+ i = -(yp n+1 / 'x^^^d^a^, 

7«A = t^ /2 c qn c lnl q~ 1/2 <£ n r n +i = {yp n+ J Xp n+1 ) 1/2 \ n c In , 

of the gauge transformed Lax parameters in terms of the points p n and the set of 2N points q„ and r n . Note that we can 
identify: 

r„ = {VvJxpJ 1 ' 2 = v^yp/xp) 1 ' 2 , (B.16) 

so that the off-diagonal elements of the T2-Lax operator are independent from the points p n and the parameter tp is related 
to our spectral parameter by: 

%ocA" 2 . (B.17) 
It is a simple exercise to derive from the above parametrization (1B.151 > the relations: 

x q n _ 3/2 Ab„ x Iri 1/2 Xd n r n+1 

— ~1 o 7 — — 1 n , (ti. IS) 

Up n Pri'n x v n +i P n 

IK = ^/a*^ Jb ! _ = _ q L/2_te2_ t SqnSin= _^_ (B . 19) 

B.3 Generalized dilogarithm functions 

Let us fix 2N + 1 points (p,qi, ...,qpj,ri, - - . ,ri\i ) £ (C 3 )' 21 ^" 1 " 1 - 1 , then we can define the following N functions: 

Yf\z n \z' n ,z' n+1 ) ee N^pJ^)^^^^!)^^/^;)^^^^/^!), (B.20) 
of the three variables (z n , z' n , z' n+1 ) £ § 3 , where: 

m f m ( s ¥«y^~ x ^ \ K/P m ( , \ _ ( ~ & ^ K+l/P m7n 

q - FjzJ " [4My§ - O ) ' T ^ (Zn+l) ~ \ s p psUa p n+1 xi - x p -) j (a - ZL) 

andz' n = q 2h ™ £ S p , while the points (p 1; ...,p N ) £ fC 3 ) N are defined bv ( IB. 13b where the (d, ctn) £ C N are fixed by 
the conditions: p p p 

< < -tF-T '^^fl'^P = l Vn£{l,...,N}. (B.22) 

Note that the Y < ^ l \z n \z / n , z' n+1 ) are well defined functions of the three variables (z n , z' n , z' n+l ) £ § 3 . Indeed, they satisfy 
the cyclicity conditions w.r.t. z' n and z' n+l thanks to the normalization functions defined in ( IB.2U while the cyclicity w.r.t. 
z n is given by the conditions ( IB. 221 ). 
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B.3.1 Recursion relations 

The generalized dilogarithm functions satisfy the recursions: 



as it simply follows from the recursions (1B.21 . Moreover, we can use the recursions (IB.7b - (IB.8b to derive 



A \ Z n<l\ z m z n+ll _ ( z n+l\ ,_ _ S/jS» ^ z n + i m 9zn 



Y™ fz \z' z' ) \ z' ■'P^+i'W" i H n --1 i n 1 jfcn „-l 



and 



where: 



_ A \ z n<l\ z m z n+l) I Z n+1 \ f _ x P n H z ™ 

yn I \ I I \ - I I ) JPnPn+lin'n i 3/q„ 

' qX^ n \nt ^n+lJ V *n / 1 — 1 z^—q- 



•/p„P„ + iq„l'n -, Wo. _ *' ' (B.25) 



Wq„pX0) ^r„p„ +1 (z(Q) 

w qn v n (z(o))w tn? (z(Q)Y 



f t = 7 (b.26) 



B.4 Baxter operator construction by gauge transformation 

Here, we present explicitly the construction of the Baxter Q-operator which can be seen as a "generalized chiral Potts" 
transfer matrix and the computation of the coefficients of the corresponding Baxter equation. 

B.4.1 Q-operator 

Let us define the kernel of the Q-operator by the ansatz: 

N 

Qa(z,z') = (z|Q A |z'> = n^'^K.Ci), (B.27) 

n=l 

where (z| and |z) are the generic elements of the left and right u„-eigenbasis and the Y^\z n \z' nl z' n+l ) are defined in 
( IB. 23b . So the recursion (IB. 23b . by using the parameter identifications ( lB.18l )-( |B.191 l, reads: 

o+l fn+1 ^ 

^\z nq \z' n ,z' n+1 ) r n+1 z' n+1 c„a„ ^ _ (l + ^7^^ 

and then: 

<z|L„(A) 21 Q A |z'> = f[ Yf (^tX^+OMA^Y^^K^^^O, (B.29) 

h^n,h— 1 

i.e. the defining condition of the Q-operator. 
B.4. 2 Baxter equation 

Let us derive the Baxter equation, from the condition (|B.29t , we have that: 



(B.28) 



(z|t 2 (A)Q a |z') = l[L n (X) 11 y[ n \z n \ Z ' n ,z' n+1 ) + l[L n (X) 22 y[ n) (z n \ Z ' n7 z' n+1 ), (B.30) 

n=l 
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where: 



UWnY^K,^) = [(Xa n - ^V)v„ - X + ^^)y- 1 ]Y x n) (z n \z' n , z' n+1 ) 

1 r n Z n r n z n 



and: 



/3„C„ A f -i , Aa„ Zn \ 



L„(A)22Y^ n) (z n |z;,z; +1 ) = [(7„/A+ 1/2 a " Z ", — )v n - (<5«A - I — ^ !L )v n 1 ]Y^ ) (z n |z;,z; +1 ) 

Acl„r z n+l 'i 



V ^™Wfl^ A+ T' ) 7: „l/2. rAQ ~, \ V n Y A ^"K^n+lJ- {B.iZ) 



A I 1 - 



'/ 



Now the recursions (lR24l - (lR25t . by (lB18t - jBJ9t , read: 

yWf, r ' ^ \V2 1 + 2^ 1+ *a 

' A l z ™9l z n' _ / z n+l \ r T /3 re r„ g 1 / 2 r„+i€„ 



Y™ (z Iz' Z' -1 V Z' y Jp " P "+i q ° r ° 1 | q^hn Zn 1i Aa n ' (D.jj, 

and 

Y A \ z n1\ z ni z n+l) _ I z n+l \ ,_ _ _ a„ z„ g 1 / 2 ^ , R ^ 

Y£>>„K,<+i) UW i.av^ ^ A d , +1 < ±1 - ^ ■ ' 

So using ( IB. 331 ) we get: 

/ , n 1/2 , i 1 _L <? 1/2A bn 

L„(A)iiY^(z„|z n ,z„ +1 ) = - ( —J- ) (3 n (q -^rX+-)— ^7 — 



/3n<E„ A 1 + 



X /p„p„ + 1 q„r„Y>/^(z n |z„, z n+1 ), (B.35) 

and analogously using ( IB. 341 ) we get: 

Z n \ a i "nlDn x , 1 n q 1 / 2 ^ 



, i x 1/2 , 1 1 - Ad " r "+ 1 

L„(A) 22 Yi")(z„|z;,z; i+1 ) = --^ (-£-) ft.fa^A+l)— 

' n+l \^7»+l/ Pn<n-n A \ — 2 2 ™ 

X /pT>Pr, + iqr, 1 "™^9A( z n| Z 7i: z 'n+l)- (B.36) 

Finally, we have the Baxter equation: 

r 2 (A)Q A = a B {X)Q x/q + d B (X)Q qX , (B.37) 



with coefficients which read: 



! ™ Hn^n 1 i T79 

Tl— 1 9 ' '"n.+l'En. 

N 1 K 1 - Ad " r " +1 

**(a) = (-i) N n ft»w^(l + ^ A )r ^fer - (BJ9) 



Remark 7. It is worth pointing out that these proofs hold also for the case of the r 2 -model on the chP curves just imposing 
that the r n in ( IB. 16b are all equal to the r in ( 15.41 ) plus the requirement that q n and r„ are on the curve. 



B.5 Connection between generalized Baxter Q-operator and SOV construction 

In the next two subsections we point out the connection among the averages of the coefficients of the Baxter equation (IB. 371 ) 
and the averages values of the Yang-Baxter operators in this way establishing the connection with the coefficients of the 
SOV-representations. 
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B.5.1 Averages of Baxter equation coefficients as eigenvalues of At (A) 

In this subsection we show that the averages of the coefficients of the Baxter equation and the cyclicity parameters 
{a\, ...,cr^} are completely characterized in terms of the matrix Ad (A) composed by the averages of the Yang-Baxter 
generators. In particular, the averages of the Baxter equation coefficients ( IB. 381 ) and ( IB. 39b coincide with the eigenvalues 
of At (A) while the eigenstates of M(A) fix the {erf, cr^}. 

\B.5.1\ a Matrix characterization of cyclicity conditions 

Let us recall that in appendix IB. 31 we have introduced the cyclicity conditions (IB. 221 ) to assure that the generalized Y- 
functions are well defined. These cyclicity conditions plus the closure condition erj^ +1 = erf fix the parameters {erf, cr^}. 
In particular, the following matrix characterization holds: 

Lemma 11. Let us define the 2x2 complex matrices: 

k n = A n ---A 2 A 1 , (B.40) 

with 

(„p „p T v q ,v _ oiP o,p t p„.p _ „p „p t p t p \ 
„p q p ~.p T p _ „.p„.p T p„.p _ „p „p T p T p j ' 
q n r n yp r n a<pBq n -^pUp °q n °r n u 'q n u 'r n / 

then 

^+i=/A„K) (B.42) 
and a\ is a solution of the quadratic fix-point equation: 

°l = fM\ CB.43) 

Here, we have defined: 

_ ax + b (a b\ 

Ja{x) = : for any 2 x 2 complex matrix A = I . (B.44) 

cx + a \cdl 

Proof. Let us point out that the cyclicity conditions ( IB. 221 ) can be seen as recursion relations on the a? n parameters. In 
particular, the cyclicity conditions for the couple (ofj, cr^ +1 ) is clearly equivalent to: 



IT 



P 

ri+1 



IaSO- (B.45) 



It is worth to recall now that (IB.44t defines a group morphism between GL(2, C) and the group of Mobius transformation: 

}ab = I a ° fs VA,BeGL(2,C). (B.46) 
Then, from this property we derive (IB .421 . which for n =N and by the closure condition cr^ +1 = a\ gives (IB. 43b . □ 

\B.5.1\ b Solution of cyclicity conditions and averages of Baxter equation coefficients 
Proposition 6. Let us denote: 



e+\ „ _ fe±\ fet \ „ , . , / e. 



At (A) \^±j = r»+ (A) , At (A) J = SI- (A) ^zj , (B.47) 
then we have one of the following cases: 

o\ = ~(x p /y p r/ 2 e+/e+, ]J a B (\q n ) = Q+ (A) , d B (\q n ) = (A) , (B.48) 

ri— 1 n—1 

of = -{x p /y p y' 2 e±/eZ, f[ a B (Xq n ) = Q. (A) , f[ d s (A<f) = Q + (A) . (B.49) 
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Proof. Let x and y be two complex numbers and define: 

\ I = A ■ \ - I (B.50) 



x' \ , I x 



\y 

it holds: 

Sa ( - ) = - r (B.51) 
Let (/ii , fci) a couple of complex numbers such that <j\ — h\jk\, then we can define the following sequence: 

(B.52) 



£ -M£. 



which satisfies the property u? n — h n /k n thanks to ( IB .42b . ( IB. 50b and ( IB. 521 ). So the closure relation (IB.431 I is equivalent to 
the solution of the spectral problem for the 2x2 complex matrix An : 

^ ft) -*ft) -fc)- 

From the definition (IB.52I I. it is simple to show that the sequence (/i n , k n ) enjoys the interesting property: 

U rrPuP — tP h 1 rr p i/P — tP 

Kn+1 = ( x P y P - x p v p ) ~E_ ft " = f x P y P _ x p w p ) " +li/ P IT" (B 54) 

from which we derive: 
being: 

p W (r p 7/P — rP i/P } a p iiP — r p 

l[a B (\ q h ) = [] fp p % V (B.56) 

P ^ (tPiiP — tP i/P ) rr p 7/P — tP 

To write the above averages formulae for the coefficients as (A) and (is (A), we have used the formulae ( IB. 381 ) and dB.39K 
the correspondence (IB. 15) and the fact that from the definition ( IB. 261) of p q r , we have: 

p— 1 / p p p p p p \ 21 

Now remarking that: 

^1 ±L= (t^ L -) = A ^ and detA N =A AN A; N , (B.59) 
with A^ n the second eigenvalue of An, we obtain: 

P N , . v -p P iV / . v -p 

nMV)=A AN n(^) , n^(^)= A A,n(^) • ^) 

h— 1 n— 1 v 7 /i— 1 n— 1 v 7 

Finally, to derive our results we have just to remark that the following identities hold: 
where C n (A) is the average matrix (13.121 ) and so by Proposition Q] 

/ N ■ \P 



m = 



H*M P/2 %(A) -Wf P ' 2 °Y (B.62) 



□ 
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B.5.2 Relations to the averages of the SOV coefficients 
The following identities: 

p p 

l[ a B (\q n ) + J] d B (Xq n ) = A(A) + P(A), (B.63) 

P P 

11 a B (Xq n )d B {Xq n ) = detM(A) = J|det q M N (Ag l ), (B.64) 



are simply consequences of Proposition [6] They are important as they explicitly imply that the sum and the product of 

the averages of the Baxter equation coefficients are Laurent polynomials in A. Let us define the following set of complex 
numbers: 

zb^U^Yu^^)} (B.65) 
where the rja were defined in ( 13.21 . then we can prove: 
Lemma 12. It is possible to fix: 

af| A£zB ^0, (B.66) 

then, under this choice, the following identities holds: 

p p 

A(A) = JJ a B (Xq n ), V(A) = l[d B (\q n ), VA e z B . (B.67) 



n=l 



Proof. Let us remark that from Proposition |6]it holds: 



v _ I , W2 *4(A) - V(A) + ey/(A(A) - V(A)) 2 - 4g(A)C(A) 

— \ x p/yp) 2 : 

with e = ±, then in the zeros of B(A), we have: 



(B.68) 



(7 



'1, 

Then it is clear that with the choice e = + it holds: 

p 



p =~{x p /y p y' 2 {A{A)-V{A)) or a{ _ = 0. (B.69) 



Ua B (\q" 



A{A), VA 6 z B , (B.70) 



and so our statement ( |B.67t follows once we use the quantum determinant relation ( IB. 64) at the zeros of B(A): 

A(A)V(A) =det q A4 N (A), VA g z B . (B.71) 

□ 

The Lemma [T2l and the characterization of subsection 13. 5 1 of the coefficients of the Baxter equation ( 13.481 implies that we 
can always chose: 

f[ A(Xq n ) = f[ a B (Xq n ), f[ u(Xq n ) = f[ d B (Xq n ). (B.72) 

n—1 ri—l n—1 n—1 

The above average identities imply that there exist two functions g (A) and / (A) such that: 

MV = J I rr-M^ d B (X) = jj^m, CB.73) 
here we show that in fact these relations are stronger being: 

9 (A) oc / (A) . (B.74) 
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The above statement follows by remarking that: 

a,B(\)dB(\/q) 



= N B (X), (B.75) 



a(A)d(A/<?) 
where: 

i \ \ "d'n^n / 1 H T72 - 1 - T72 / 

and from the following: 

Lemma 13. The function Nb (A) is the identity for all the representations of the Ti-model. 

Proof. Let us rewrite the functions /p„f>„ +1 q n r„ by using their definition (IB. 26b and the parametrization (IB. 18t -( |B. 19b ; then 
they reads: 



_ l a n+ i r n a n p n \ ^ t g g 0n - - „ ? ,g„---«-n 

and we can write: 



n B (x) = n 



1 l + 9 -a*fll/2Sii_i_ l-g2fc g -l/22» Ar ' 

-1 y " € n an 



(B.78) 



1 + ^ /2 (fr^) P i-9 p/2 (frAr„) P 

Thanks to the cyclicity conditions ( IB. 221 ) the r. h. s. of the last equation is 1 as announced. □ 

Remark 8. Let us point out that for general representations the generalized Baxter Q-operator is not proven to define 
a commuting family w.r.t. the spectral parameter A as well as it is not proven to commute with T2(A). In this general 
setting the eigenstates of T2(A) are not necessarily eigenstates of the generalized Baxter Q-operator. Vice versa, the Baxter 
equation (IB. 37b implies that any eigenstate of the generalized Baxter Q-operator is also eigenstate of T2(A) which so will 
be still constructed as explained in Theorem [3] Finally, let us remark that the previous results on the coefficients implies 
that the Baxter equations characterizing the generalized Baxter Q-operator and the Baxter Q-operator constructed by SOV 
coincide up to a gauge transformation. 

C. Properties of the cofactors Qj(A) 

Lemma 14. Let t(X) be a real Laurent polynomial of degree N in A, then the cofactors of the matrix D{\) satisfy the 
following identities: 

c h+l , k+l {\) = c M (A<f) Vi,Me{i,...,2l + i}, (C-i) 

Ci,i(A)=Ci,i(-A), C 2 ,i(A) =<7 N C!, 2 (-A). (C.2) 

and: 

(Ci,i(A))* = Ci,i(eA*), (Ci, 2 (A))* = Ci j3I +i(eA*). (C.3) 

Proof. The proof of properties (IC.lt and (IC.2) coincides step by step with that given in Lemma 4 of 111. Let us prove the 
property (1C.31 for the cofactor Ci.i (A) = det2/-Di,i(A), where: 

£>i,i(A) = \\t(Xq h )S Kk - z(\q h )6 h , k+1 - d(Ag h )<y h , fc _i|| 1 < h < 2Jil < it < 2J , (C.4) 

then by the properties under complex conjugation of t(A) and being a(A)* = e N d(eA*) it holds: 

(Di,i(A))* = ||e N (t(e\*q h )6 htk - d(e\*<T h )5 h , k+1 - a(AV^)<^-i) \\ 1<h<2l 1<k<2l > (C5) 
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Let Di i be the 21 x 21 matrix of columns: 



Ca 1 ' 1 = Cfi; 1 , Vae{l,...,2l}, (C.6) 



where denotes the column a of the matrix X, and let D l \ be the 21 x 21 matrix of rows: 



where denotes for the row a of the matrix X, then the following identity holds 



Ra l11 = R°l£, Vae {!,..., 2/}, (C.7) 



det^; 1 K (eA*) = (det J D 1 , 1 (A))* (C Xl i(A))* s Ci,i(eA*). (C.8) 



Analogously, it holds: 

N 



(det£> (1 , 2) , (1 , 2) (A))* = e N det Z) (1 , 2) , (1 , 2) (eA7g), (C.9) 



and by using the expansions: 

2/ 



Ci, 2 ,+i(A) - JJa(Ag fc ) + d(A/g)detX> (lia)l(1 , 2) (A/ 9 ) 1 (CIO) 
Ci, 2 (A) = JJd(A« fc ) + a(A ? )det£) (li2)i(lj2) (A). (C.ll) 

h=l 

the second identity in dC.3l > is just a consequence of iC.9\ . □ 
Lemma 15. For any t(X) 6 £ r2 the following identities are verified: 

Zci,i(A) n Z Cl , 2 (A) =^Ci,i(A) n 2 Cli21+1 (A), (C12) 

one/ 

3 s pM c Z Cl l(A) ^ 5p, Ao nZ Cl 2(A) ^ 0, (C.13) 
forS p ^x = (A ,gA , ...,q 2/ A ) anyp-string. 

Proof. The proof of identity (IC. 12b follows step by step the proof given in Lemma 5 of the paper (TJ. Let us assume that 
dC. 1 3I > is not satisfied, then ( 14.171 ) implies that S p \ C Z Cl 2(+1 (a) from which ( |4.20t holds only if S p Ao C Z d( - A ), which is 
not verified by the definition of d(A). □ 
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